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Substituted oligothiophenes have a long history in the field of 
organic electronics, as they often combine outstanding elec- 
tro-optical properties with the ease of synthesis. To assist the 
rational selection of the most promising structures to be syn- 
thesized, there is the demand for tools that allow prediction of 
the properties of the materials. In this study, we present strat- 
egies for synthesis and computational characterization, with re- 
spect to the fluorescence behavior of oligothiophene-based 
materials for organoelectronic applications. In a combined ap- 



proach, sophisticated computational methodologies are direct- 
ly compared to experimental results. The M06-2X functional in 
combination with the polarizable continuum model in a state- 
specific formulation for excited-state solvation proved to be 
particularly reliable. In addition, a semiclassical approach for 
describing the vibrational broadening of the spectra is em- 
ployed. As a result, a robust procedure for the prediction of 
the fluorescence spectra of oligothiophene derivatives is pre- 
sented. 



1. Introduction 

Organic electronics have gained much interest from industry in 
recent years, and this field covers the use of organic semicon- 
ductors within electronic devices, for example, OFETs (organic 
field-effect transistors), OLEDs (organic light-emitting diodes), 
and OPVs (organic photovoltaics). 111 The application of these 
"plastic electronics" in the field of display technology enables 
the construction of very-thin, flexible, energy-efficient, light- 
weight panels that are highly suitable for customer electronic 
displays. 121 Considering the benefits of this technology, very in- 
tense academic and industrial research activities are on the 
way to overcome deficits of these materials (e.g. limited life- 
times). Thus, the quest for novel materials is an ongoing pro- 
cess, especially because OLEDs are believed to have enormous 
economic potential in the near future. 

However, the synthesis of novel organic compounds can be 
an expensive procedure, considering both the costs of fine 
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chemicals, reagents, and catalysts and man power. Thus, it 
would be a great benefit to obtain information on the physical 
properties of interest of the substance prior to synthesis. On 
the basis of such knowledge, the most promising candidates 
can be selected within a substance class under investigation. 
In the case of OLED materials, luminescent properties are of 
particular interest. As the molecules used are, in general, of 
considerable size, the availability of reliable and computation- 
ally efficient methods for the prediction of photoluminescent 
properties is essential. Poly- and oligothiophene-based com- 
pounds are widely used within the field of organic electron- 
ics. 131 They display remarkable conductivity with excellent 
chemo- and thermostability, and they are, in general, readily 
synthesized or structurally modified to meet electronic and 
spectroscopic requirements. Hence, one of our ongoing re- 
search topics involves the modification of the substance class 
BRA-nT (Figure 1). 

A small subset, namely, BMA-nT, was studied by Noda et al.™ 
and found to be suitable for the fabrication of single- and mul- 
tilayer OLEDs. Moreover, it was shown that the use of triaryl- 




a) 



b) 



Figure 1 . a) Substance class BRA-nT (BMA-nT: R = CH 3 , n = 1 -4; studied by 
Noda et al.). b) Model compounds DMA-nT (n = 1-3). 
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amine subunits had a positive influence on both the conduc- 
tivity and the tendency to form amorphous glasses. Thus, com- 
binations of (oligo)thiophenes with triaryl moieties are widely 
used in all fields of organic electronics. 151 

For BRA-nT, variation of both the conjugation length n of 
the oligothiophene linker as well as the R substituents of the 
triarylamine are interesting, especially with regard to the ef- 
fects of the applied modifications on the electrooptical charac- 
teristics. Emission wavelengths were calculated by time-depen- 
dent density functional theory (TDDFT), 161 as recent develop- 
ments in this area allow the reliable and efficient treatment of 
molecules of considerable size, which makes it the method of 
choice for this study. Solvent effects were taken into account 
by applying a polarizable continuum model (PCM). Aside from 
its standard linear response formalism for excited states, the 
state-specific solvation model was also used. 171 Moreover, to 
understand and describe spectral broadening, which is an im- 
portant factor in, for example, the design of white OLEDs, we 
included vibrational information. Simulations based on the 
Wigner distribution of the vibrational ground state of the first 
excited electronic state proved to be appropriate to perform 
this task. As a result of this project, a reliable procedure can be 
presented for the calculation of the entire fluorescence spectra 
of oligothiophene-based compounds. 

2. Results and Discussion 

As already mentioned, bis-substituted derivatives containing 
(oligo)thiophene substructures have gained considerable at- 
tention in the field of organic functional materials [e.g. BHA- 
1T ra (4); R = phenyl, n = 1; Scheme 1]. As a consequence, 




3a-c: n = 1-3, R = -CH 3 
4:n = 1,R=™(3 



Scheme 1. Synthesis of 3a-c (DMA-nT) and 4 (BHA-1T). Reagents and condi- 
tions: i) 1 a-c (1.0 mmol), boronic ester 2a,b (3.0 mmol), degassed /PrOH/ 
H 2 0, KOfBu (3.0 mmol), (IPr)Pd(aHyl)CI (2.0 mol% for 3 a-c, 5.0 mol% for 4), 
reflux («80°C), 1-2 h. 



a large number of oligothiophene derivatives have been syn- 
thesized and characterized. Additionally, these types of com- 
pounds have been subjected to theoretical studies. These 
cover, for example, investigations of vibrational (IR and 
Raman) 191 and electrochemical 1101 properties as well as the treat- 
ment of excited states. 1111 However, in most cases no experi- 
mental data was supplied for comparison with the calculated 



properties. In contrast, our study focuses on the application of 
modern and computationally affordable methods to molecules 
relevant in the field of materials science, as we combine our 
theoretical results with experimental data of well-characterized 
compounds. As a future perspective, these data can be inter- 
esting for subsequent investigations, for example, to provide 
reference data for other approaches. To minimize the computa- 
tional effort during the initial stage of the project, simplifica- 
tions were applied to BRA-nT. Replacement of the relatively 
large triarylamine system with a considerably smaller dimeth- 
yl(4-aminophenyl) group leads to model compounds 3a- 
c (DMA-nT, n=1-3; Scheme 1), which were used as mimics to 
benchmark the applied computational techniques and parame- 
ters used within the calculations. Compound DMA-1T was pre- 
viously characterized in 1998 when obtained as a byproduct 
during the synthesis of 1,2-dithiin derivatives. 1121 However, a reli- 
able synthetic protocol has not been published so far. A selec- 
tive approach toward DMA-nT and BHA-1T could be developed 
through Suzuki cross-coupling strategies by utilizing dibromo- 
thiophenes 1 a-c and dioxaborolanes 2a,b (Scheme 1). 

The best results were achieved by using the N-heterocyclic 
carbene (NHC) catalyst (IPr)Pd(allyl)CI, [131 which afforded the de- 
sired materials in good to excellent yields. According to our ex- 
perimental experience, this catalytic system, in general, shows 
high reactivity for the benzene-thiophene linkage. Crystalliza- 
tion of BHA-1T from cyclohexane yielded single-crystalline 
samples; the molecular structure of 4 is illustrated in Figure 2. 




Figure 2. Molecular structure of 4 (BHA-1T). 

Relative to other excited-state methods, the TDDFT approach 
is favorable because of its low scaling behavior with respect to 
increasing size of the system, which leads to acceptable com- 
putational costs, even for large systems. However, the results 
obtained can strongly depend on the functional used. 1141 Stan- 
dard functionals may particularly suffer from overstabilization 
of charge-transfer excited states 1151 and are not reliable for 
large conjugated systems. 1161 To select an appropriate function- 
al for the description of the fluorescence properties of the sub- 
stance class under investigation, three different functionals (i.e. 
B3LYP, 1171 PBE0, [6a ' 181 M06-2X [6bl ) were considered. With these 
functionals, benchmark calculations were performed on the 
model DMA-nT compounds and solution-based computational 
results were compared to the experimental values. 

In Figure 3, the optimized geometries of both the ground 
state (S 0 ) and the first excited state (S,) of DMA-1T are present- 
ed. A significant structural change is observed, and the most 
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Figure 3. Optimized geometries of DMA-IT (3a) (M06-2X/SVP, gas phase); 
left: S 0 , right: S,. 

obvious difference is the torsion angle between the phenyl 
substituent and the linking thiophene unit. Taking the M06-2X- 
optimized structures as an example, this angle is about 24° in 
S 0 , whereas it decreases to 0 G in S,, which thus leads to 
a planar excited-state geometry. By comparing the optimized 
geometries of S 0 and S„ the same planarization can be ob- 
served for DMA-2T and DMA-3T (for illustrations see the Sup- 
porting Information, Figures S1 -S3). When analyzing the bond 
length alternation in DMA-1T, one can recognize a significant 
contraction of the C 1 C 1 ' thiophene bond as well as the C2— 
C3 linkage, whereas the C1— C2 bond is elongated. Moreover, 
the quinoidal distortion of the benzene rings is enhanced. 
Both phenomena, the planarization and the quinoidal distor- 
tion, are processes that are well known for systems of this 
t yp e neb,i9] Q eta j| ec | information on the bond lengths in both 
the S 0 and S, geometries of DMA-1T is presented in Table 1. 



Table 1. 


Bond lengths of DMA-IT (M06-2X/SVP, gas phase). 


Bond 


Length [A] 








So 


s, 




ci-cr 


1.4222 


1.3788 




C1-C2 


1.3747 


1 .4234 




C2-S 


1.7356 


1 .7640 




C2-C3 


1.4678 


1.4162 




C3-C4 


1.4018 


1 .4270 


4 ''~\\ 


C4-C5 


1.3877 


1.3786 




C5-C6 


1.4123 


1.4192 




C6-N 


1.3767 


1.3738 




C7-N 


1.4418 


1.4424 





In a first step, gas-phase calculations for the vertical emis- 
sion from the minimum geometry of the first excited state (S,) 
were performed (Table 2). First, it may be noticed that by in- 
creasing the number of repetitive thiophene units (i.e. going 
from DMA-1T to DMA-3T) and, therefore, elongating the it 
electron system, the fluorescence wavelength is also increased. 
By examining the results for the vertical emissions obtained by 
the three different functionals under investigation, decreasing 
wavelengths are found moving from B3LYP to PBEO to M06-2X, 
which is similar to the trend reported in ref. [20]. For the small- 
est system (i.e. DMA-1T), this difference between the excitation 
energies obtained from B3LYP and M06-2X calculations 
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Table 2. Calculated 


gas-phase 


vertical emissions at various 


levels of 


theory. All data denoted in nm. 








Level of theory 


Vertical emission 








DMA-IT 


DMA-2T 


DMA-3T 


BHA-1T 


B3LYP/SVP 


440(461) 


508(533) 


561(590) 


- 


(6-31 1++G**) 










PBE0/SVP 


428(446) 


490(513) 


539(564) 




(6-31 1++G**) 










M06-2X/SVP 


407(426) 


457(480) 


492(516) 


431(450) 


(6-31 1++G**) 










RI-CC2/TZVPP 


409(416) 


455 


487 


448 


(QZVPP) 











amounts to a shift of 33 nm (0.23 eV). In the case of DMA-3T, 
this difference is larger at 69 nm (0.31 eV). A comparison with 
the RI-CC2/TZVPP reference (Table 2) and the experimental re- 
sults (see below) clearly favors the M06-2X results. When com- 
paring the excited-state geometries obtained from the three 
different functionals, a nearly ideal flat structure is derived 
from B3LYP and PBEO computations, whereas M06-2X leads to 
a slightly curved geometry (Figure 4). 




Figure 4. Optimized S, geometries of DMA-IT (3 a, top) and DMA-3T (3 c, 
bottom); in both graphs: upper (bent) structure calculated with M06-2X/SVP, 
flat structure with B3LYP. 

All these data suggest that B3LYP (and to some extent also 
PBEO) is insufficient for the description of the conjugation ef- 
fects in extended jt systems. The insufficiency of standard 
functionals in such cases has been reported before, 1161 and it 
has been pointed out that inclusion of exact Hartree-Fock (HF) 
exchange is decisive in this context. 1251 Thus, it can be under- 
stood that the high nonlocality functional M06-2X provides the 
most stable results. In a similar sense, it was reported that the 
range-separated CAM-B3LYP functional outperformed the PBEO 
functional in the case of delocalized excited states. 1201 To esti- 
mate basis set effects, the TDDFT calculations were carried out 
with the enhanced 6-311 ++G** basis set. These calculations 
led to a general modest redshift of about 25 nm (0.1 eV) for all 
systems and functionals considered. However, considering that 
this shift is rather small and uniform, it appears that the SVP 
basis set is sufficient, given that the resulting excitation ener- 
gies are generally somewhat higher than those in the com- 
plete basis set limit. DMA-1T was also treated at the RI-CC2/ 
QZVPP level, and the results showed a slight decrease in exci- 
tation energy relative to the TZVPP results. In summary, at this 
stage it can be concluded that the M06-2X functional is in ex- 
cellent agreement with the RI-CC2 reference. To test the applic- 
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ability of the developed strategy on "real-life materials" and to 
obtain a rough approximation of the computational time, we 
included one derivative of substance class BRA-nT (Figure 1), 
namely, BHA-1T (4; R = H, n=1), in our performance study to 
serve as a representative example. In this case only, the M06- 
2X functional was considered and good agreement with the 
RI-CC2/TZVPP reference was obtained. Comparison to DMA-1T 
shows that the outer phenyl groups have a small effect on the 
fluorescence, and they result in a redshift of about 25 nm 
(0.17 eV). 

Selected molecular orbital plots at the M06-2X level are 
shown in Figures 5 and 6; a larger set of orbitals, also including 



LUMO 



HOMO 



Figure 5. Molecular orbital plots (DMA-IT, S„ M06-2X/LR-PCM, THF). 




LUMO 



HOMO 



Figure 6. Molecular orbital plots (DMA-3T, S„ M06-2X/LR-PCM, THF). 



information on the other functionals, is given in Section C of 
the Supporting Information. In all cases, the main contribution 
to the excited state is a transition between the highest occu- 
pied molecular orbital (HOMO) and the lowest unoccupied mo- 
lecular orbital (LUMO). These are always distributed over the 
whole planar system, which also includes the phenyl rings and 
the terminal amino nitrogen atoms. In the cases of larger sys- 
tems, contributions of the outer phenyl rings to the LUMO are 



somewhat diminished. The HOMO of BHA-1T also contains 
some contributions from the outer phenyl rings. When com- 
paring the shapes of the orbitals between the three functionals 
for the DMA-nT series, it is observed that these are virtually in- 
distinguishable. However, an interesting difference between 
the functionals is observed when considering the weight of 
the HOMO/LUMO transition. In calculations with the B3LYP and 
PBEO functionals, the HOMO/LUMO transition contributes with 
more than 99% for the whole DMA-nT series (for B3LYP, this 
value was in most cases even above 99.9%). By contrast, in the 
case of M06-2X this value steadily decreased along the series 
giving 98, 97, and 95% for n=1, 2, and 3, respectively, with 
the HOMO-1/LUMO + 1 transition coming into play. By com- 
parison, the RI-CC2/TZVPP values are 95, 95, and 94%, respec- 
tively. Some of us have recently discussed such properties, 1211 
which are related to the concept of natural transition orbi- 
tals. 1221 In particular, it should be noted that additional transi- 
tions coming into play are related to electron correlation 1231 
and entropy. 1241 An interpretation in the present context is that 
a pure HOMO/LUMO transition represents a homogeneous ex- 
cited state, or in other words, an unbound exciton. Additional 
structures arising from electron-hole interactions can only be 
included through additional independent transitions. In this 
sense, it is expected that for bound excitons, the weight of the 
HOMO/LUMO transition should decrease as the size of the 
system is increased. Specifically, it has been pointed out that 
the amount of HF exchange is important for providing bound 
excitons. 1251 In agreement with this interpretation, the best re- 
sults are obtained with M06-2X containing, at 54%, a signifi- 
cantly larger proportion of HF exchange relative to that of the 
other functionals. In summary, we interpret our results in the 
sense that for the system sizes considered here (i.e. up to five 
conjugated rings), new excitonic effects start to come into 
play. In spite of their good performance for smaller molecules, 
the B3LYP and PBEO functionals are incapable of describing 
these correctly. By contrast, M06-2X appears to provide a very 
good description for molecules of the size considered here. 
Furthermore, for even larger conjugated systems, inclusion of 
full HF exchange (at least after range separation) may be nec- 
essary. 

As the next step, the fluorescence spectra in solution were 
considered. These data are, in particular, interesting, as it has 
been shown (e.g. for BMA-nT) [4] that the photoluminescence 
spectra of solutions are similar to those recorded for the fabri- 
cated OLED. To investigate the influence of the solvent, three 
different solvents, namely, n-heptane (HEP), dichloromethane 
(DCM), and tetrahydrofuran (THF) were selected. The fluores- 
cence spectra in these solvents were measured. Computation- 
ally, solvent effects were included at the PCM level in its linear 
response (LR) and state specific (SS) formulations (see Compu- 
tational Details). The main effect, a redshift resulting from an 
increase in the solvent polarity, could be reproduced with all 
computational models. 

Further results can be extracted from Table 3:1) The above- 
mentioned trend in emission wavelengths (B3LYP> PBE0> 
M06-2X) remains true if solvent corrections are applied; to il- 
lustrate this, the LR-PCM vertical emissions calculated by using 
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Table 3. Comparison of calculated vertical emissions by applying SS-PCM 
(LR-PCM results in parentheses) in different solvents and experimental 
(exp.) values. All data denoted in nm. 

Parameter 



Solvent 



HEP 



DCM 



THF 



Vertical emission [a] 


Functional 


A IT 

UrvIA- 1 1 


UIVIA-z 1 


UMAo 1 


BHA-1T 


B3LYP 


452(46/) 


C~i A i C A A\ 

Dz4p44) 


COT 1 £Lf\A\ 
Z>OZ(OU4) 




r tstu 


43o(4o:S) 








M06-2X 


41 7(429) 


470(486) 


506(523) 




exp. 


413-437 


n.a. Ibl 


n.a. Ibl 




B3LYP 


471 (505) 


555(596) 


624(667) 




PBEO 


458(489) 


536(572) 


599(635) 




M06-2X 


434(461) 


495(526) 


534(567) 




exp. 


438-458 


490-511 


525-548 




B3LYP 


470(502) 


552(592) 


621(663) 




PBEO 


457(486) 


534(569) 


595(631) 




M06-2X 


433(457) 


493(523) 


532(564) 


455(474) 


exp. 


429-451 


481-507 


517-543 


444-468 



[a] All experimental spectra showed two distinct maxima, [b] Substance 
insoluble in n-heptane. 



the three different functionals are superimposed on the experi- 
mental spectrum in Figure 7 (top). 2) M06-2X, in particular with 
SS-PCM, provides very reliable results; values obtained from 
both B3LYP and PBEO are redshifted by several tens of nano- 
meters for both LR-PCM and SS-PCM corrections. 3) For all 
functionals, LR-PCM overestimates the solvatochromic redshift, 
whereas the SS-PCM results are very close to the experimental 
values. A more-detailed discussion of this overestimation and 
additional information on these solvation models are given in 
ref. [19c]. Furthermore, it is worth mentioning that oscillator 
strengths (f) increase quite strongly when moving from DMA- 
IT to DMA-3T (e.g. M06-2X/SS-PCM in THF: 1.40, 1.77, and 2.13 
for DMA-1T, DMA-2T, and DMA-3T, respectively). In all cases, 
the transition dipole vector lies in the molecular plane and 
points along the aromatic backbone. 

For a graphical illustration of the obtained accuracy, Figure 7 
(bottom) shows the experimental fluorescence spectrum of 
DMA-1T (THF solution) as well as the M06-2X/SVP-based theo- 
retical values obtained when applying no (gas phase) LR-PCM 
or SS-PCM for the influence of the solvent. The gas-phase 
value is strongly blueshifted relative to both maxima of the ex- 
perimental spectrum. It lies close to the onset of the spectrum, 
that is, the experimental value for the 0-0 transition, which 
should lie at significantly higher energies than the vertical fluo- 
rescence. By contrast, the LR-PCM value is redshifted to such 
an extent that it is situated at higher wavelengths with respect 
to both experimental emission maxima. The most accurate pre- 
diction for the vertical emission is obtained when the SS-PCM 
correction is applied. Considering Table 3, this method yields 
a value that is situated between the two maxima of the experi- 
mental spectrum in all cases except DMA-1T in DCM, in which 
it is just below them. Figures 8 (top) and 9 show the experi- 
mental spectra and the calculated vertical emissions for DMA- 
nT and BHA-1T, respectively, highlighting the suitability of the 
applied computational models. 



DMA-1T (Exp., 370nm Exc.) 
DMA-1T (M06-2X, THF) 
DMA-1T (PBEO, THF) 
DMA-1T (B3LYP, THF) 




500 550 
Wavelength /nm 




— DMA-1T (Exp., 370nm Exc.) 

— DMA-1T (M06-2X, g.p.) 
•-■ DMA-1T (M06-2X, THF: LR-PCM) - 
•-- DMA-1T (M06-2X, THF: SS-PCM) 



500 



550 



Wavelength /nm 

Figure 7. Vertical emissions of DMA-IT in THF calculated by using the B3LYP, 
PBEO, and M06-2X functionals and LR-PCM as solvent correction (top). Verti- 
cal emissions of DMA-IT calculated by using the M06-2X functional in the 
gas phase and in THF solution by applying LR-PCM or SS-PCM (bottom). Ex- 
perimental spectra measured from THF solutions. 



The next step in the course of this study was to examine the 
band broadening of the fluorescence spectrum. In addition to 
recording the experimental spectra, simulations were per- 
formed by a semiclassical approach as outlined in the Compu- 
tational Details Section. For reasons of computational efficiency 
and ease of implementation, spectra simulations were per- 
formed by applying a LR-PCM. As outlined above, this leads to 
a slight redshift in the simulated spectra. However, the band 
shape should be affected to a lesser extent. The experimental 
spectra for the DMA-nT series are illustrated in Figure 8 (top). 
All three spectra are similar in appearance with a slightly no- 
ticeable vibronic structure featuring peaks that are, in all cases, 
separated by about 1400 cm" 1 . This feature probably arises 
from aromatic stretching modes, considering that the corre- 
sponding internal coordinates are strongly altered during the 
excited-state structural relaxation (cf. Table 1). Such modes 
around 1400 cm" 1 were indeed found in the ground-state 
normal mode analysis of these molecules. In contrast to the 
fluorescence spectra, the S, peaks in the absorption spectra of 
these molecules did not show any vibronic structures (see Sec- 
tion D of the Supporting Information). This difference most 
probably originates from the planarization in the excited state, 
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T 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■ 

350 400 450 500 550 600 650 700 750 800 
Wavelength /nm 



T 1 1 1 1 1 1 ' 1 1 1 1 1 1 r 




1 1 ■ 1 ■ 1 1 1 1 [ 1 1 1 1 1 1 1 

350 400 450 500 550 600 650 700 750 800 

Wavelength /nm 



Figure 8. Experimental (top) and simulated (M06-2X, LR-PCM, bottom) emis- 
sion spectra of DMA-nT (n = 1-3) in THF. Vertical emissions (M06-2X, SS-PCM) 
are indicated by dashed lines. 

as described above, and a concomitant increase in the torsion- 
al barrier (cf. ref. [19b]). The full widths at half maximum 
(FWHM) for the experimental spectra are 0.40, 0.38, and 
0.39 eV for DMA-1T, DMA-2T, and DMA-3T, respectively. The si- 
mulated spectra (Figure 8, bottom) are in acceptable agree- 
ment with the experimental ones with respect to the position 
of the maximum as well as the broadening of the band. How- 
ever, the FWHM broadening of these spectra (0.55, 0.53, and 
0.48 eV) is somewhat too large relative to the experimental 
ones. A notable difference can be seen at the blue edge, at 
which point the experimental spectra decrease more sharply 
than their computational counterparts. These differences may 
be related to the treatment of torsion, which is, by default, 
simply considered as a harmonic oscillator in Cartesian coordi- 
nates. Considering that, in our semiclassical approach, only the 
modes of the planar excited-state structure are included explic- 
itly, this still gives reasonable results with the agreement 
shown above. By contrast, in a Franck-Condon approach, in 
which the wavefunctions of low-frequency torsional modes 
over multiple nonplanar wells in the ground state also have to 
be computed, the situation becomes more difficult and re- 
quires special techniques or artificial planarization of the mole- 
cu l e [i9c,26] ^ r jg 0rous approach in which all of the torsions are 
considered as hindered rotors in curvilinear coordinates is, in 
principle, also possible in our approach. However, such a treat- 



ment was already quite involved for the smaller bithiophene 
system, which contains only one torsion, 1190 - 261 and therefore 
this treatment is beyond the scope of this work. 

It should be noted that the semiclassical approach used 
here is only capable of reproducing the spectral envelope 
without vibronic progressions. The structure that is visible in 
these spectra is probably only of stochastic nature and would 
disappear with increased sampling (requiring significantly en- 
hanced computational effort). However, the presented compu- 
tational protocol carries the potential to provide interesting 
additional information in the event that either decomposition 
of the spectra is performed (cf. ref. [27]) or the involved modes 
are analyzed in more detail. In particular, it would be of great 
interest to understand which modes comprise the principle 
components of spectral broadening, considering that such in- 
formation could be used to modify the spectral properties in 
a rational design process. Such an analysis, which is planned as 
a future project, should also shed new light on the differences 
with regard to the experimental results at the spectral edges. 

Figure 9 summarizes the whole protocol as applied to BHA- 
1T. In this case, good agreement is again found. The vertical 
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Figure 9. Experimental and simulated (M06-2X, LR-PCM) emission spectra of 
BHA-1T in THF. Vertical emission (M06-2X, SS-PCM) is indicated by a dashed 
line. 



emission (M06-2X, SS-PCM) is located between the two 
maxima of the experimental spectrum, which is reproduced by 
the spectra simulation over a wide range, except for the blue 
onset (M06-2X, LR-PCM). 

3. Conclusions 

The aim of this study was to establish a computational proto- 
col for the prediction of the electronic emission spectra of 
functional organic molecules. Oligothiophene-based com- 
pounds, widely applied in the field of organic electronics, were 
characterized with regard to their fluorescent properties. Three 
representative model compounds were synthesized, analyzed, 
and used to benchmark some critical parameters for the calcu- 
lations. Beyond various functionals under investigation, M06- 
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2X gave highly accurate results in comparison to both experi- 
ment and RI-CC2 benchmark calculations, whereas for B3LYP 
and to a lesser extent PBEO, the values were redshifted, espe- 
cially when applied to molecules with extended conjugated jt 
systems. A tentative explanation of this phenomenon was 
given in terms of excitonic structure starting to arise in sys- 
tems of this size. Solvent effects were included by applying 
a PCM in its linear response and state specific formulations, 
and the latter provided superior results. An accuracy of ap- 
proximately ± 15 nm can be expected for the prediction of ex- 
perimental fluorescence maxima in solution by calculating ver- 
tical emissions. The same satisfactory results can be obtained 
also for larger molecules of immediate interest, for example, 
BHA-1T. Moreover, vibrational broadening was simulated in so- 
lution by utilizing the framework of the Newton-X software. 
The obtained spectra are generally in good agreement with ex- 
perimental data. However, an overestimation of the blue onset 
could be observed, and this in turn leads to a somewhat in- 
creased spectral width. Investigation into possible reasons 
leading to this effect will be a matter for future studies. 

Having verified the utility of this computational protocol, 
the next task will be to obtain more-detailed insight into the 
underlying physics, and its influence on the spectral properties 
of these compounds, which will be achieved by extended anal- 
ysis of our data. Aside from considering influences on the posi- 
tion of the absorption maximum, a particular focus will also be 
on the molecular modes affecting spectral width. On the basis 
of this understanding, new molecules with desired properties 
can be designed and synthesized and their properties deter- 
mined. We believe that the close connection between synthe- 
sis, spectroscopy, and computation, which is provided by our 
laboratory, provides an ideal environment for this feedback 
cycle. 



Computational and Experimental Details 

General Methods 

Density functional theory (DFT) and time-dependent (TD) DFT cal- 
culations were performed by using three different functionals: 
B3LYP,' 171 PBE0, [6a ' 181 and M06-2X. [6bl Most TDDFT calculations were 
performed with the polarized double-^ SVP basis set. 1281 In the 
cases specified, the larger 6-311 ++G** basis set 1291 was used as 
a reference. Geometry optimizations were performed without sym- 
metry constraints and all ground- and excited-state minima were 
verified by normal mode analysis. Solvent effects were included 
through the polarizable continuum model (PCM) 1301 in its linear re- 
sponse (LR-PCM) 1311 and state specific (SS-PCM) 171 formulations, 
always considering the equilibrium time regime (eq) for the excited 
state. Vertical fluorescence energies were obtained after geometry 
optimizations in the S, state at the TDDFT level, both in the gas 
phase and by using LR-PCM (eq). For benchmark purposes, higher- 
level calculations with the use of the second-order approximate 
coupled cluster method (CC2) [321 with resolution of the identity ap- 
proximation 1331 (RI-CC2) were performed as well by using the 
TZVPP and QZVPP basis sets. 1341 Subsequently, vibrational broaden- 
ing of the spectra was simulated by using a Wigner distribution 
formalism, as described in ref. [35], by applying a phenomenologi- 
cal broadening of (3 = 0.1 eV. Taking the output of the frequency 



computations as the source of the normal modes, the Newton-X 
framework is capable of fulfilling this task very effectively. Knowl- 
edge upon the excited-state normal modes is used to generate 
a set of nonequilibrium geometries with a selectable number of 
data points, which in these particular cases was 500 (DMA-nT) or 
250 (BHA-1T). At each point, the S,^S 0 excitation energy was cal- 
culated by using the M06-2X functional at the TDDFT level of 
theory. Considering that these are all independent single-point 
energy calculations, this step can be easily performed in parallel by 
distributing them to a large number of compute nodes. After com- 
putation of all emission energies, the data is collected to construct 
a vibrationally broadened spectrum. In the case of solution-phase 
spectra, not only the emission energies but also the geometries 
and frequencies were computed in solution by applying LR-PCM 
(eq). In addition to the semiclassical spectra, a Franck-Condon (FC) 
simulation of the emission spectrum of DMA-IT was attempted in 
Gaussian 09. This led to the problem of an almost vanishing over- 
lap integral between the two vibrational ground states of the S 0 
and S, states, and thus the calculation was aborted. A solution to 
this problem requires either an artificial planarization of the rele- 
vant torsional coordinates or a more-involved treatment of these 
coordinates (cf. [ ' 9c ' 2a ) and was not attempted. All TDDFT computa- 
tions were performed using the Gaussian 09 package, revision 
A.02. [3a RI-CC2 calculations were carried out with Turbomole 6.3. 1371 
Spectra simulations were performed in the framework of Newton- 
X [35,38] ^ n j nter f ace 0 f Newton-X to Gaussian 09, which was devel- 
oped in the course of the project, is available in the current 
Newton-X release. 

If not stated otherwise, chemicals were obtained from commercial 
sources and used without further purification. Column chromatog- 
raphy was performed on silica 60 (Merck, 40-63 |im) by using dis- 
tilled solvents as given. The following substances were synthesized 
according to the literature: 5,5'-dibromo-2,2'-bithiophene (1 b), [391 
5,5"-dibromo-2,2':5',2"-terthiophene (1 c), [40] N,N-dimethyl-4-(4,4,5,5- 
tetramethyl-1,3,2-dioxaborolan-2-yl)benzeneamine (2a), [411 and N,N- 
diphenyl-4-(4,4,5,5-tetramethyl-1,3,2-dioxaborolan-2 L)benzenea- 
mine (2b). [41M21 

NMR spectra were recorded at 400 MHz for 'H and at 100 MHz for 
13 C with a Bruker Avance DRX-400 spectrometer. Data for 'H NMR 
are reported as follows: chemical shift in parts per million from tet- 
ramethylsilane with the residual solvent signal as an internal refer- 
ence (CD 2 CI 2 : 3 = 5.32 ppm, [D,]TFA: (3 = 11.50 ppm; TFA = trifluoro- 
acetic acid), multiplicity (s = singlet, d = doublet, t = triplet, and 
m = multiplet), coupling constant, and integration. "C NMR spectra 
are reported in ppm from tetramethylsilane by using the central 
peak of the solvent as reference (CD 2 CI 2 : 5 = 54.0 ppm, [D,]TFA: 
5 = 116.6 ppm), multiplicity with respect to proton (deduced from 
APT experiments, s = quaternary C, d = CH, t = CH 2 , q = CH 3 ). Laser 
desorption ionization (LDI) high-resolution (HR) mass spectra were 
obtained with a Synapt HDMS (Waters, U.K.). Experimental fluores- 
cence characterizations were performed with an Edinburgh FLS920 
fluorometer; signal processing was realized by using Savitzky- 
Golay methods. Single-crystal X-ray diffraction data were recorded 
with a Bruker Kappa APEXII CCD diffractometer, A(MoKa) = 
0.71073 A. Refinement of F 2 was done by using the SHELXTL pro- 
gram (version 2008, Bruker AXS, Inc., Madison, Wl). 

Syntheses 

General procedure for the synthesis of DMA-nT and BHA-1T com- 
pounds: KOfBu was added (3.0 mmol) to a suspension of dibromo- 
thiophene 1 a-c (1.0 mmol) and boronic ester 2a,b (3.0 mmol) in 
degassed /PrOH/H 2 0 (v/v = 3/1, 12 mL). Subsequently, under an 
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argon atmosphere a solution of (IPr)Pd(allyl)CI [13al (2.0mol% for 
3a-c, 5.0 mol% for 4) was injected into the stirred mixture by sy- 
ringe and heated to reflux (f»80 °C). Upon completion of the reac- 
tion (monitored by TLC, 1-2 h), the products were purified by chro- 
matography or crystallization. Analytical samples of 3a-c and 4 
were (re)crystallized from pyridine or cyclohexane, respectively, 
prior photophysical characterization. 

4,4'-(2,5-Thiophenediyl)bis[N,W-dimethylbenzeneamine] (3a, DMA- 
IT): The reaction mixture was poured into H 2 0/CHCI 3 (1:1), and the 
aqueous layer was extracted with CHCI 3 . The combined organic 
layers were dried with Na 2 S0 4 and concentrated under reduced 
pressure. Flash chromatography (cyclohexane/DCM, 6^100%) on 
triethylamine-preconditioned silica gel yielded 3 a (241 mg, 75%) 
as an orange-brown solid. 'H NMR (400 MHz, [D,]TFA): (3 = 7.98 (d, 
J = 8A Hz, 4H), 7.68 (d, J = 8.4Hz, 4H), 7.52 (s, 2H), 3.53 ppm (s, 
12H); 13 C (100 MHz, [D,]TFA): r5 = 150.5 (s), 142.8 (s), 126.7 (d), 123.4 
(s), 122.0 (d), 113.0 (d), 40.8 ppm (q); HRMS (LDI-qTOF): calcd for 
C 20 H 22 N 2 S [M] + 322.1504; found 322.1632. 

4,4'-[2,2'-Bithiophene-5,5'-diyl]bis[N,W-dimethylbenzeneamine] (3 b, 
DMA-2T): The reaction was performed following the general proce- 
dure. Upon completion, the reaction mixture was filtered through 
a sintered funnel, and the residue was washed with /PrOH/H 2 0 
(3:1), EtOH, and light petroleum ether and then dried under 
vacuum. Crystallization from pyridine yielded 3 b as a light orange 
solid (295 mg, 73%). 'H NMR (400 MHz, CD 2 CI 2 + [D,]TFA): 6 = 7.79 
(d, J = 8.6 Hz, 4H), 7.58 (d, J = 8.6 Hz, 4H), 7.37 (d, J = 3.9 Hz, 2H), 
7.28 (d, J = 3.9Hz, 2H), 3.33 ppm (s, 12H); ,3 C (100 MHz, CD 2 CI 2 + 
[D,]TFA): (3 = 141.3 (s), 141.1 (s), 138.4 (s), 137.1 (s), 128.0 (d), 126.4 
(d), 126.0 (d), 121.5 (d), 47.8 ppm (q); HRMS (LDI-qTOF): calcd for 
C 24 H 24 N 2 S 2 [M] + 404.1381; found 404.1393. 

4,4'-[2,2':5',2"-Terthiophene-5,5"-diyl]bis[A/,A/-dimethylbenzenea- 
mine] (3 c, DMA-3T): The reaction was carried out on a 0.5 mmol 
scale; the workup was performed in a manner identical to that for 
3 b. Crystallization from pyridine afforded 3 c as an orange solid 
(160 mg, 66%). 'H NMR (400 MHz, CD 2 CI 2 + [D,]TFA): 5 = 7.78 (d, J = 
8.6 Hz, 4H), 7.57 (d, J = 8.6 Hz, 4H), 7.36 (d, J = 3.9 Hz, 2H), 7.24 (d, 
J = 3.9Hz, 2H), 7.21 (s, 2H), 3.33 ppm (s, 12H); "C (100 MHz, 
CD^ + IDJTFA): (5 = 141.2 (s), 140.8 (s), 138.6 (s), 137.1 (s), 136.7 
(s), 128.0 (d), 126.4 (d), 125.6 (d, 2C), 121.5 (d), 47.8 ppm (q); HRMS 
(LDI-qTOF): calcd for C 28 H 26 N 2 S 3 [M] + 486.1258; found 486.1165. 

4,4'-[2,5-Thiophenediyl]bis[N,A/-diphenylbenzeneamine] (4, BHA-1T): 
Workup was performed in analogy to that performed for 3 a. Flash 
chromatography (cyclohexane/DCM, 5^10%) with silica gel yield- 
ed 4 (526 mg, 92%) as a bright yellow solid; single crystals were 
obtained from cyclohexane. 'H NMR (400 MHz, CD 2 CI 2 ): 6 = 7.50 (d, 
J = 8.4Hz, 4H), 7.32-7.24 (m, 8H), 7.20 (s, 2H), 7.15-7.01 ppm (m, 
16H); 13 C (100 MHz, CD 2 CI 2 ): (5 = 148.0 (s), 147.8 (s), 143.1 (s), 129.9 
(d), 128.9 (s), 126.7 (d), 125.1 (d), 124.2 (d), 123.7 (d), 123.7 ppm (d); 
HRMS (LDI-qTOF): calcd for C 40 H 30 N 2 S [M] + 570.2130; found 
570.1812. Crystal data (yellow prism of 0.40x0.30x0.20 mm from 
cyclohexane): C 40 H 30 N 2 S, FW = 570.72, 7"= 296(2) K, triclinic, space 
group Pi, a = 9.8178(2) A, 6 = 12.3368(2) A, c= 1 2.9148(2) A, a = 
88.365(1)", fi = 76.765(1)°, y = 78.962(1)°, V= 1494.34(5) A 3 , Z=2, 
Pcakd = 1 -268 gcnrr 3 , ^ = 0.141 mm -1 . Of 31978 reflections measured 
up to 6 = 30°, 8659 were unique. Refinement of F 1 concluded with 
R1 =0.0659 and w/?2 = 0.1450 for 388 parameters and all data. 
CCDC-858839 contains the supplementary crystallographic data for 
this paper. These data can be obtained free of charge from The 
Cambridge Crystallographic Data Centre via www.ccdc.cam.ac.uk/ 
data_request/cif. 
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